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BY HYBRID COMPUTERS 

By Raymond S. Lim and William D. Cameron 
Ames Research Center 


SUMMARY 


This paper describes a method for analyzing power and cross -power spectra 
from 0.1 Hz to 20 kHz with general purpose hybrid computers. The method is 
valid for analyzing continuous and discrete spectra of periodic, aperiodic, 
complex, and stationary random signals. A spectrum is analyzed entirely on 
the analog computer by a pseudo band -pass filter method. The bandwidth and 
the averaging time are varied to suit the computation requirement of each 
spectrum. The digital computer is used for automatic control of the analog 
computer, for data acquisition, and for compensation of magnetic tape static 
skew. The definition of power and cross-power spectra, and the validity of 
various methods of analysis are comprehensively reviewed. 


INTRODUCTION 


Power and cross -power spectra are analyzed by filtering, phase shifting, 
squaring, and averaging. During the past 20 years, many methods have been 
devised for investigating problems in subjects ranging from aeronautics to 
zoology and many special analyzers have been built. Reference 1 presents a 
survey of spectrum analyzers. The flexibility of such analyzers is limited in 
frequency range, filter bandwidth, signal skew error compensation, averaging 
time, and the number of channels that can be analyzed simultaneously. 

Recently, digital computers have been used for spectrum analysis (in a very 
limited range) . 

The expanding technology of aeronautical and biomedical engineering 
requires a more sophisticated method for spectrum analysis. This requirement 
prompted the development of hybrid computers for spectrum analysis. The pur- 
pose of this paper is (l) to present a hybrid -computer technique for analyzing 
continuous and discrete spectra of periodic, aperiodic, complex, and station- 
ary random signals from 0.1 Hz to 20 kHz, (2) to clarify well-established 
theories that describe random data, and (3) to clarify mechanization methods 
for power and cross-power spectrum analyses. 

Existing hybrid computers, with a few additional special instruments, can 
be programmed to form a high quality spectrum analyzer. Hybrid computers are 
useful for spectrum analysis because they have the speed and flexibility of an 
analog computer, the dynamic range and accuracy of a digital computer, and the 
automatic output of results for documentation and plotting. 
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At Ames Research Center, a hybrid computing system is used to compute 
power and cross -power spectra of all types of data. The data, recorded on 
analog magnetic tape, have a frequency range from 0.01 Hz to 20 kHz, and may 
or may not he stationary throughout the entire data record. The length of the 
data record ranges from 50 msec to 60 sec with a dynamic range to 100 dB. 

Three channels of data can he analyzed simultaneously hy the analog computer. 
The digital computer automatically controls the analog computer, performs the 
data acquisition, and compensates for the static skew error of the magnetic 
tape. It should be noted that any random data recorded on magnetic tape can 
he treated as stationary random data if the entire record length is analyzed. 


NOTATION 


Ci 2 (m) 

cospectrum 


QisO 0 ) 

quadspectrum 


Rll(r) 

autocorrelation function 


Ri2(t) 

cross -correlation function 


T 

averaging time 


X(w) 

Fourier transform of x(t) 



x(t) passed through a hand- 

-pass filter with bandwidth 

x L (t) 

x^)(t) delayed hy 90 ° 


x(t) 

time average of x(t) 


x(t) 

ensemble average of x(t) 


Aj 

bandwidth 


$n(w) 

power spectral density of 

Xl(t) 

$22( w ) 

power spectral density of 

x 2 (t) 

**» 

power spectral density of 

x(t) 

w c 

cutoff frequency of low-pass filter 

0) 

o 

local oscillator frequency 
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ANALOG SYMBOLS 





amplifier 



integrator 



track/store amplifier 


P in, 


P\r\ 2 



+ ^in, P in 2 
K(cu) 


K(o >) = 100 at dc 


multiplier 


P in 


{T) k(u>V in 


potentiometer 



" A *," 1 


X wc (I) 


band -pass filter 


low-pass filter 


X (t) 


AVE X(t) 


average circuit 




MATHEMATICAL THEORY 




Definition of Power Spectrum 


There are at least two ways to calculate the power spectrum of stationary 
random data: (l) by direct calculation and (2) by transforming the autocorre- 

lation function. By these methods, three basic mathematical definitions for 
the power spectrum are generally derived, only two of which are valid for 
random data. Let $n(w) be defined as the power spectral density (psd) ; 
x ( t ) , a stationary random signal; and T, the interval of the data to be 
analyzed. The three definitions of psd are then: 


1. From direct calculation, 

3>n(w) = lim 

T->oo 

2. From direct calculation. 


$n(u)) = lim — 
T— >co X 


x(t)e ^ dt 


“T 


x(t)e dt 


-T 


3- From transformation of the autocorrelation function, 


^n(w) = / Rn(T)e dr 


= 2/ Rh(t)cos wt dT 
J o 


2^T / x L (t)dt 

Aw^o ,J o 


- lim 


where 


Rii(t) = lim - / 

T -*00 Jq 


x(t)x(t - T)dr 


( 1 ) 


( 2 ) 


(3a) 

(3b) 

(3c) 


(4) 


Equations ( 3 ) are mathematically correct and are excellent definitions 
for mechanizing on either digital or analog computers. For digital computation, 
equation (3b) is generally used, while for analog computation, equation (3c) 
is generally preferred. Equation (3a) is mathematically related to the corre- 
lation analysis of random data. The proof for equation (3a) is sometimes 
called Wiener's theorem for autocorrelation, and it can he found in refer- 
ence 2. The derivation of equation (3c) from (3a) is relatively 
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* straightforward, but since this derivation is not often found in published 
literature, it is developed in appendix A. Equation ( 3 c) can be implemented 
to obtain the density spectrum as follows: 

1. x(t) is passed through a band-pass filter with bandwidth Au) to 
obtain x^(t), 

2. x A Jt) is squared through a multiplier to obtain x AuA), 

3. finally, x Aw (t) is averaged to obtain 

Equations (3) are applicable to stationary random signals as well as periodic 
and nonperiodic signals. Equation (3c) is an estimate of equation ( 3 a). 

Reference 3 (Aseltine) shows that equation (2) is correct for x(t) . 

(The wavy bar denotes the ensemble average.) This equation has little 
practical use in spectrum analysis. 

The first definition, equation (l), has been shown to be mathematically 
invalid for all classes of random signals (ref. 3, 4, 5, or 6). 


Definition of Cross -Power Spectrum 

Cross -power spectrum analysis obtains the amplitude and phase -spectrum 
information between two independent random signals, Xi(t) and x 2 (t) . From the 
electrical engineering point of view, a cross -power spectrum analyzer is a 
phase meter for measuring the phase and amplitude relationship between two 
complex or random signals. The cross-power spectrum is a vector quantity, or 
a complex function. The real part is called the cospectrum and the imaginary 
part , the quadspectrum. 

The cross-power spectrum function ^12(^0) can be defined as the Fourier 
transformation of the cross -correlation function in a manner similar to defin- 
ing the power spectrum as 


$12 ( w ) 



Ri 2 (t) 


e _JwT 


dT 


(5) 


Mathematically, this is the only known correct definition for random signals. 
The definition of the cross-power spectrum function from direct calculation as 




r T . , 

lim jj ! xi(t)e^ W dt 
■p— >00 Uq 



( 6 ) 


is incorrect when xi(t) and X2(t) are random signals (for the reason given in 
the previous section for power spectrum). However, if xi(t) and X2(t) are 
periodic functions, Papoulis (ref. 7) shows that equation ( 6 ) is correct. 
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Equation (5) can be defined in terms of the cospectrum and the 
quadspectrum as 


where 


and 


^12( U ) = Ci 2 (w) + <3Qi2(w) 


Qi 2 (w) = 


-r 

Ri 2 (t) 

-00 


= lim 
Ad )-»0 

1 

2AwT ^ 

T->oo 



R 12 ( T ) 

-00 


= lim 

1 

2AwT 


>T 


X W t)x2 Aw (t)dt 


sin ut 




Xl A 3 t)x2 Aw (t)dt 




( 7 ) 

(8a) 

(8b) 


( 9 a) 


( 9 b) 


In the above equations, Ci 2 (u>) is the cospectrum, Qi2(w) is the quadspectrum, 
Ri2(t) is the cross-correlation function, x 2 ^(t) is x 2 (t) passed through a 
band -pass filter with bandwidth Aw, and xg^Jt) is shifted 90°. 

(Equations (8b) and ( 9 b) are derived in appendix B.) 


Basic Method of Power Spectrum Analysis 


x (t) 

Aw 

x AqjM| 

V / 

,171 



/ \ 

\ 



<X> xx (aj)Aa> 


Figure 1.- Analog method of power spectrum 
analysis . 



Figure 2.- Analog method of cross-power 
spectrum analysis. 


The analog method of computing 
power spectra is well-known. As indi- 
cated by equations (3), (8b), and (9b), 
the basic method of computing the power 
and cross -power spectra is by band -pass 
filtering, phase shifting, squaring, 
and averaging the data (figs. 1 and 2, 
respectively) . 

Actually, the methods in figures 1 
and 2 cannot be reasonably mechanized 
because it is difficult to design 
high-frequency narrow band-pass fil- 
ters (Ah of 1 Hz at 10 kHz) , band- 
pass filters that can automatically 
select the bandwidth and the center 
frequencies independently, and two 
band -pass filters with an accurately 
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matched phase for cross -spectrum analysis. Consequently, the basic analog 
* method is not generally used. In its place, a pseudo band -pass filter (some- 
times called tracking filter) is mechanized by the heterodyne principle. But 
because of this filtering technique, the method is valid only for analyzing 
data with a sufficiently long record. It is inaccurate for analyzing short 
records of a transient nature (such as impact research or heartbeat data), 
unless the data are made periodic. (This will be clarified in subsequent 
discussions.) In the present technology, analog multipliers (for multiplying 
and squaring) and 90 ° phase -shift networks with the required bandwidth are 
readily available. 


TYPES OF RANDOM SIGNAL TO BE ANALYZED 


The types of random data to be analyzed, by the techniques described in 
this report, can be classified into three groups: 

1. Aerodynamic data on nonsteady phenomena. The frequency range of 
interest is from 10 Hz to 20 kHz with a dynamic range of 7^ dB and a data 
record 60 sec long. The data are similar to white noise mixed with a few 
periodic signals, and may or may not be stationary within the 60 -sec length. 
The spectrum is continuous. 

2. Impact research data. The frequency range of interest is from 1 to 
2500 Hz with a dynamic range of 80 dB and a data record 50 to 100 msec long. 
The data are similar to a decaying oscillation mixed with a small amplitude 
random signal. The spectrum is not continuous. These types of data must be 
converted to periodic data before the analysis so that they may be averaged. 

3 . Physiological and biological data. The frequency range is from 0.1 
to 1000 Hz with a dynamic range of 60 dB and a data record between 0.1 to 

2 sec long. Most physiological data analyzed are heartbeats and brain waves, 
both human and animal. The spectrum is not continuous. These data must also 
be converted to periodic data before the analysis. 


HYBRID COMPUTER METHOD OF SPECTRUM ANALYSIS 


Method of Power and Cross-Power Spectrum Analysis 

Power and cross -power spectra are analyzed by a hybrid computer in a 
manner similar to the method discussed under Basic Method of Power Spectrum 
Analysis. There are many combinations of analog and digital computers that 
are called hybrid computers. The hybrid computer used to implement the spec- 
trum analysis discussed in this report is the EAI HYDAC 2000 system, which 
consists of a 231R-V analog computer and a DOS -350 logic computer. The method 
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of analysis is shown in figure 3, and 
is thoroughly analyzed in appendix C. 
As shown in appendix C, the spectra of 
xi(t) and x 2 (t) are 


2 

^n(w) = 2^ (Output 1 ) 

(10a) 

2 

$> 22 (w) = 2^ (Output 2 ) 

(10b) 


Equations ( 10 ) are not ambiguous and 
are correct forms for random signals as 
well as complex or periodic signals. 

The division in equations ( 10 ) is 

for bandwidth normalization, and should be used for a continuous spectrum only. 
For cross-power spectrum, the cospectrum and quadspectrum are 

Ci 2 (w) = (Output 3 ) (lla) 


Figure 3.- Method of power and cross-power 
spectrum analysis. 


Qi 2 (w) = - 2 ^* (Output k) 


In particular, for inputs of 

xi(t) = Ai sin ojit 
x 2 (t) = Ai sin( Ul t - 0 2 ) 

0 2 = variable, from 0 ° to 360°, Ci 2 (io) and ^12(10) are 


Ci 2 (oj) - 


Qi 2 (oj) - 


A x 2 B 2 


c°s 0 2 1 ^ 


Ai 2 B 2 


Sin 


(lib) 


(12a) 

(12b) 


regardless of 101 > oj 0 or uji < aj Q . This means that there is no ambiguity 
on the output sign of the quadspectrum. 

The implementation of the power and cross -power spectrum analysis between 
computer elements of the method shown in figure 3 I s 


1 . Analog computer; all spectrum computations. 

2 . Digital computer; 

(a) Logic computer - for automatic control of the analog computer 
and analyzer output digitizing and storage. 
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(b) IBM 70^0/709^ - off-line spectrum scaling, coherence function 
computation, cross -power phase -angle computation, magnetic tape recorder 
static-skew error compensation, and plotting. 

In addition to the standard HYDAC 2000 system, a programmable two -phase 
oscillator system and a low-pass filter system are also required. 

Perhaps the best way to explain how the hybrid computer is used to 
mechanize the spectrum analysis is to divide the system in figure 3 into six 
major parts, and then describe each part in detail; these six parts are: 

(l) a programmable two-phase oscillator, ( 2 ) a pseudo band-pass filter, 

( 3 ) squaring and averaging, ( 4 ) output equations, ( 5 ) automatic program con- 
trol and analog data digitizing, and ( 6 ) scaling and static skew-error 
correction. 



Figure 4.- Programmable two-phase oscillator. 



Figure 5-- Auxiliary equipment required for 
spectrum analysis. 


Programmable two -phase oscillator . - 
A block diagram of the programmable 
two -phase oscillator circuit is shown 
in figure 4. The hardware is shown in 
figure 5 • The frequency of the spec- 
trum to be analyzed is punched onto 
the paper tape to be read by the block 
reader. The frequency synthesizer is 
a new instrument and is simply a digi- 
tal programmable oscillator. The two- 
phase shifter is an all-pass active 
network with two outputs. One output 
is the reference phase; the other is 
the -90° phase. The two output phases 
are always 90 ± 5 ° apart for the fre- 
quency within the band of 0.1 to 
5000 Hz. Before 1966, a two-pnase 
synthesizer was not available. The 
amplitude response of the best single - 
phase frequency synthesizer is ±12 per- 
cent, and this larger amplitude 
variation must be compensated in the 
two -phase shifter to within ±1 percent 
or better. At the present, it is rea- 
sonable to expect that a two -phase 
synthesizer with an amplitude response 
of ±2 percent or better and a phase 
linearity of ± 2 ° throughout the fre- 
quency range will soon be available. 
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Pseudo band-pass filter . - The mechanization of the pseudo hand -pass 
filter is shown in figure 6. Consider the pseudo hand -pass filter 1, which 
consists of multiplier 1 and low-pass filter 1. The input signal x(t) is 
multiplied, or heterodyned, with the 0° phase of the oscillator. The resul- 
tant output is a sum and difference frequency pair for each frequency Wj_ of 
x(t) . At the output of low-pass filter 1, only those frequencies within w c 
are passed. Thus the multiplier low-pass filter combination is, in effect, a 
pseudo hand-pass filter. The center frequency is determined hy the oscillator, 
and the bandwidth is equal to 2w c . In effect, is x(t) passed 

through a hand -pass filter with bandwidth Ajo. 


The synthesized hand-pass filter in figure 6 is not a true hand-pass 
filter. It measures only that component of the signal inphase with the refer 
ence oscillator, as shown in appendix C. Therefore, the application of one 
pseudo hand -pass filter for spectrum analysis is not sufficient. Two pseudo 
hand-pass filters must he used to measure both the inphase and out -of -phase 
components of the signal, since, in general, there is a random phase angle 


PSEUDO BAND-PASS FILTER NO. I 



PSEUDO BAND-PASS FILTER NO.2 

Figure 6.- Pseudo band-pass filter and ripple 
cancellation. 


between the input data and the oscil- 
lator. If this phase angle is repre- 
sented hy 9 ± , the output of the first 
pseudo hand-pass filter is 

cos(ojt + ) and the second, 

Aj_ sin(wt + Q±) . To obtain a true mea- 
sure of the signal, it is necessary to 
add these two components vectorially in 
the manner of 


A-j 2 = [ Aj_ cos(urt + 0q)] 2 

+ [ Aj_ sin(u)t + @i) ] 2 


PSEUDO BAND-PASS 

3O0K|| FILTER 


PREFILTER 



Figure 7*" Analog mechanization of pseudo 
band-pass filtering system. 


The pseudo hand-pass filtering 
circuit is shown in figure 7- Only the 
Xi(t) channel is shown; the X2(t) 
channel is the same. The circuit shown 
in figure 7 should he self -explanatory. 
The outputs of low-pass filters 1 and 2 
are xi^(t) and x W t) ' respectively. 

The phase and the amplitude responses 
of the xi(t) channel and the Xa(t) 
channel must he matched in order for 
the cross -power spectrum outputs to he 
meaningful. Any mismatch will he 
system error. The system error in 
using the HYDAC 2000 system and the 
auxiliary equipment is approximately 
±4° of measurement with a full scale 
of 36O 0 . 
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Squaring and averaging . - The squaring and averaging circuit system is 
shown in figure 8 for $n(.w), and in figure 9 for Qi 2 (w) . The circuits for 
$ 22 (m) and Ci 2 (w) axe assumed to he understood. True integration is used. 

The integration time is set and controlled by the DOS- 35 O logic computer in 
accordance with the mode -logic table of the integrator (A05) • A05 is pro- 

grammed to integrate for a fixed length of time T and then to hold its value 
at the end of time T. During this hold period, track/store amplifier A06 
tracks A05 and stores its value. After AO 6 has stored the value of A05, A05 
is commanded to reset. Operation of the circuit in figure 9 is similar to 
that in figure 8 . 




Figure 8 . - Squaring and averaging system, 
Xi(t) power spectrum channel. 


Figure 9*“ Squaring and averaging system, 
Qi 2 ( w ) channel. 


Output equations . - There are four outputs from the analyzer when two 
channels of data are being analyzed: hi ( M ) > $ 22 (w) , Ci 2 (w) , and Qi 2 (w) . 

After the signal is traced through the analog circuits and all gain factors 
are accounted for, these four outputs become 


On((h) 


^ 22 ^) 


C J.2 ( h) ) 


Qi2(^) 


2 (Output 1) 1 S A 

(13) 

" 7.29xl0 6 KfiKi2 ^ S 0 

2( Output 2) 1 S A 

(14) 

~ 7 . 29 X 10 ^! k| 2 £** S o 

2( Output 7) 1 S A 

(15) 

“ 7 - 29 Xl 0 6 K 11 K£ 1 K 12 K 22 S o 


2 (Output 4) 1 S A 

( 16 ) 

7 . 29Xl0 6 KiiK 2 iKi2K 2 2 


where the above equations are in mean-square power , and 


magnetic tape loop speed playback when analyzing data 
magnetic tape speed at which the data were recorded originally 


11 




The factor S^/S Q is necessary because the multipliers in the 23IR-V are 
accurate only to about 5 kHz. To analyze to 20 kHz, a speed reduction of 4 is 
required. The normalization of the output equations by the bandwidth is valid 
only for continuous spectra. (A continuous spectrum is one whose output is 
directly proportional to the analyzer bandwidth.) 

Automatic program control and analog data digitizing . - The complete 
program for controlling the analog computer for automatic spectrum analysis 
is provided by the DOS-350 logic computer. The DOS-350 also digitizes the 
four outputs from the analyzer and punches these outputs on paper tape for 
subsequent spectrum scaling and static skew-error compensation on the 
IBM 70^/709^ system. Specifically, the DOS-35O performs the following 
functions during spectrum analysis: 

1. Senses the control signal recorded on one channel of the tape loop 
and starts the analysis for that frequency point as programmed on the fre- 
quency synthesizer. The splice problem on the loop is eliminated by using a 
step -type control signal on the tape loop. 

2. Generates a second set of control signals for spectrum analysis at 
higher frequencies if the data are stationary throughout the record. 

3. Generates a control signal to advance the block paper-tape reader by 
one block (one frequency point) at the end of the analysis time T. 

4. Generates integrator and track/store amplifier control signals to 
control the analog spectrum computation in a fully automatic mode. 

5. Records the number of frequency points analyzed and shuts off the 
program at the end of the analysis. 

6. Punches the program on paper tape in the proper format consistent 
with the IBM 7040/7094 system. 

Scaling and static skew-error compensation . - Scaling is defined here as 
the conversion of the spectrum outputs to engineering units of the original 
experiment. For example, if the spectrum outputs should be in psf/Hz (pounds 
per square foot per hertz) , then the original calibration of the transducer 
must be used to multiply the spectrum outputs at each frequency point. This 
type of scaling or multiplication is ideally suited to the digital computer. 

There are two ways of compensating static skew error: 

1. By an analog delay line at the output of the magnetic tape recorder, 

and 

2. By computing 0 ±2rtfT at the digital computer after scaling, where 

t is the time difference between recorder channels, and 0 is the arctangent 
of Ql2( u ) over Ci2(w) in radians. If xi(t) leads X2(t) by t, 0 - 2rtfT 
is used; otherwise, 0 + 2jtfT is used. 
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Static skew-error compensation is mandatory at high frequencies for the 
‘'cross-power spectrum since present recorders can have static skew errors as 
large as 30 (isec, which represents 216° at 20 kHz. 


Dynamic Range, Resolution, Bandwidth, Averaging Time, 
and Frequency Scan Rate 

The dynamic range of the analyzer must be carefully considered before the 
type of analog computer is selected for spectrum analysis. Dynamic range here 
is defined as: 

. . maximum spectrum output obtainable 

dynamic range = mlninwm s ^ ctr u'm output' observable 

For a high quality +100-V computer, the dynamic range is about 100 to 
0.05 or 66 dB. To obtain the maximum 66 -dB dynamic range at the output, two 
gain controls are used in each channel of the analyzer. In the xx(t) channel, 
Kxi is the prefilter gain and Kxz is the postfilter gain. Before each anal- 
ysis, a search is required for the maximum spectrum output. After that partic- 
ular frequency point is located, K n and K 12 are adjusted for a 100-V output. 
In a good ±100-V analog computer, the output is actually linear up to ±120 V 
or more. 

Resolution and bandwidth are functions of the low-pass filter setting. 

The bandwidth should be chosen in accordance with the type of data to be 
analyzed. To simplify changes of bandwidth, standard low-pass filters are 
used. 


The problem of averaging time is solved by using true averaging. The 
frequency scan rate has no meaning in this hybrid method of spectrum analysis, 
since a sweep-type oscillator is not used. After completing a frequency-point 
analysis, the next frequency point can be analyzed as soon as the system 
transient decays to zero. 


Statistical Uncertainty 

The predominant source of error or uncertainty in an analyzed spectrum is 
the random statistical variations of the sample record. Perhaps the most com- 
plete discussion, but not necessarily the easiest to understand, on the sub- 
ject of statistical error is given by Chang (ref. 8). Some error discussions 
were also given by Bendat and Piersol (ref. 9)- Confidence limits can be 
defined as relative for statistical error comparison. In reference 8, the 
confidence limits are defined as 


1 

1 + a 


< 


$(u)) 

®(w) a 


< 


1 


1 - a 


( 17 ) 


13 


vhere 


4>(w) 

®(u) a 


a 


true spectrum 
analyzed spectrum 



1 

s/SoT 

N number of sample records 

k_ a constant depending on the specified probability., tabulated as 

1 follows: 


P k p 


0.5 

0.477 

.8 

.906 

• 9 

1.163 

.95 

1.386 

• 99 

1.82 

• 999 

2.32 


Example: The signal to be analyzed has the following analyzer parameters: 

N = 1 , T = 10 sec , = 10 Hz 

Now determine the 90-percent confidence limits. 

Solution: The rms per unit error of each record is 


e = 


n/aj T 


= 0.1 


For P = 0.9 or 90 percent , kp = 1.163. Therefore, 

a = I.I63 X 0.! X = 0.165 


1 - a 


1.198 , 


1 + a 


= O.858 


Thus, there is a 90-percent probability that the true spectrum lies between 
O.858 and I.I98 times the measured spectrum. 
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Typical Results 



Figure 10.- Analysis of a sine wave. 


• MEASURED 

i. — ACTUAL 

o). . _ _ ...... 

II | *■- « 

<X> 22 (w), X* 

0R 5 

C. 2 (U0 • 

0 1 1 J - 


1 

Q| 2 (W) 

.5 


0 5 10 15 20 25 

FREQUENCY, kHz 


Figure 11.- Analysis of white noise from 
10 Hz to 20 kHz. 


a>„M 



I 250 500 


FREQUENCY, Hz 

Figure 12.- Analysis of a heartbeat pulse. 


Figure 10 shows the analysis of a 
single 0.707 Vrms sine wave at 1 kc as 
an input to both the xi(t) and x 2 (t) 
channels. The result shown in fig- 
ure 10 is correct since the analysis of 
a single sine wave for a frequency hand 
across the analyzer bandwidth is the 
amplitude response of the analyzer 
band -pass filter. 

The analysis of a band-limited 
white noise from a General Radio 
GR-I39OB generator is shown in fig- 
ure 11. The noise is first passed 
through a low -pass filter with an 
effective bandwidth of 20 kHz. The 
amplitude of the noise is then adjusted 
to read 0.45 Vrms at the output of the 
filter (i.e., the input to the ana- 
lyzer) . This gives a power spectral 
density of 10 " 5 V 2 /Hz. The scatter of 
the measured values about the actual 
values is within ±2 percent. 

Figure 12 shows an analysis of one 
heartbeat from medical research with a 
0.2-Hz analyzer bandwidth from 0.5 to 
125 Hz. 

Figures 13 and 14 show the power 
and the cross -power spectra of two 
typical 20 -kHz wind-tunnel data. These 
data are from transducers 1 inch apart , 
stations -18.19 and -17.19* These two 
figures are plotted by the digital com- 
puter side by side on a 31 X 31-inch 
sheet. The accuracy and confidence of 
these figures are confirmed by three 
CPRMS numbers. The CPRMS (wind tunnel) 
is the RMS value of the data measured 
at the wind tunnel during an experiment. 
The CPRMS (loop output) is the RMS 
value of the data after it is trans- 
ferred onto an analog loop recorder for 
spectrum analysis. The CPRMS (psd 
area) is the RMS value of the data com- 
puted in the digital computer by taking 
the square root of the area underneath 
the power spectrum curve. Under ideal 
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Figure 13 . - Power spectral of two typical 
20 kHz wind-tunnel data. 


Figure 14.- Cross-power spectrum of two 
typical 20 kHz wind-tunnel data. 
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Figure 15.- Analysis of the transfer function. Figure l6.- Analyzed transfer function of 

a human pilot. 


conditions when no error occurs, these three numbers should be equal. As 
indicated in figure 13, the error is about 3. 15 percent. 

Figure 15 shows the analysis of a known transfer function, Yp. This 
same transfer function was analyzed independently in the digital computer by 
the BOMM-Tukey method (ref. 10) . The results obtained are in close agreement 
with the actual frequency and phase response. Figure 1 6 shows the analysis of 
a human pilot. Note the close agreement of the results obtained by the three 
methods - the hybrid method, the BOMM-Fourier method, and the BOMM-Tukey 
method. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Sept. 27, 19 66 
. 124-11-04-06 
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APPENDIX A 


DERIVATION OF ® xl (u) FOR RANDOM SIGNALS 
FROM AUTOCORRECTION 


The purpose of this appendix is to show that the power spectrum 




Rii(r)e dr 


can he estimated by 


$n(w) 


= lim 

o 


1 

2&oT 


T->*oo 



(Al) 


(A2) 


To proceed with the derivation, the inverse transform of equation (Al) is 
taken. 


Rn(r) 



$ii((o)e^ WT 


doo 


(A3) 


Since 

Rn(o) > |Rh(t)| (a4) 

the rnaxirnum psd is obtained for t = 0. That is, 


Rn(o) 


0 1 i(o)) dw 


(A5a) 


= lim 

T->oo 


1 

T 


sT/2 


x 2 (t)dt 


-T/2 


(A5b) 


In the derivation of a cross spectrum, r cannot equal zero because Ris(O) is 
not necessarily the maximum -value. For this reason, the cross spectrum is 
generally complex. Equations (A5) should be clear intuitively since the total 
power of x(t) is equal to the mean square value, and also to the total area 
underneath the power spectrum curve. Now, for small dw, doj approaches Zu. 
In the limit, as doj approaches zero. 


IT 



(A6a) 


00 

= lim 2 ) ®n A (oo)Ajl) ( A6"b ) 

Ad->o ^ 
n=o 

The interpretation of equations (A6) is that integrating ®ii(w) is equivalent 
to obtaining the area underneath the spectrum curve. This area can also be 
obtained by dividing the area into n small areas with width Aj and summing 
over n. In view of equations (A6) , equation (A5a) becomes 

OQ 

Rn(0) = lim 2 ) $n. (u>)Aj (A7) 

Ao-»o L 
n=o 

or 

Rii^Co) = lim 2<Hi Aj (a))£w (A8) 

'“ aj5 AjJ-^o 


lim / $ii(w)dco = lim 
dw^o ~-oo d w-^o 


' n n 


(co) dw 


n=-oo 


from which 


= lim ^ R i^ (0) 

o 


(A9) 


or, in general notation, 

T 

®n(<*>) = lim 2^T f x L (t)dt (A10) 

A.i)->o ^ o 

>p-KX) 

where x^t) is interpreted as x(t) passed through a band-pass filter with 
bandwidth Zko. This completes the derivation. 
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APPENDIX B 


DERIVATION OF <S 12 (u)) FOR RANDOM SIGNALS 
FROM CROSSCORRELATION 

The purpose of this appendix is to show that the cross -power ’spectrum 


$12(03) - / Ri2(r)e 


-JU)T 


dT 


(Bl) 


can be estimated by 


3> 12 (w) - Cx 2 (co) + jQi 2 (w) 


(B2) 


where 


T 

C 12 (m) = lim ^ (t) Xz (t)dt 

A.r>n 


AjlT>0 

T->oo 


(B3) 


Qi 2 (oj) = lim 


T-^oo 


2AwT 




(B4) 


and xj (t) is x 2 ^(t) delayed 90 ° • To proceed with the derivation, the 

inverse transform of equation (Bl) is taken. The variable delay t cannot be 
set equal to zero because Ri 2 (o) is not equal to nor, in general, greater 
than Ri 2 (t). It is felt that the derivation that follows is heuristic, 
rather than rigorous. Now, 


r uu 


(B4a) 


[Ci 2 (w) + jQi 2 (w) ]e^ dm 


(B4b) 


Since an arbitrary function can be decomposed into a sum of an even and an 
odd function (e.g., ref. 7)> then. 

Ri 2 ( t) = Ri 2 (t) s + R12 ( t) (B 5) 
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where Ri 2 (t) s and Ri 2 (t) q are the even and odd parts of R 12 (t) , respec- 
tively. Now, the result of equating equations (B5) and (Bhb) is 


Rl 2 (T) e + R12 ( t ) q 


Cl2( (jJ ) e '^ WT &C0 


+ (j ) / Ql2(^)s^ U)T dGJ 


(B6) 


With some intuition, let us equate the even and odd parts of equation (B6) . 
That is. 


Ri 2 (r) e = 


C 12 (cj)e' :)U)T dco 


(B7a) 


Ri2(t)q = (j) 


Qi2(w)e'^ CJT dw 


(B7b) 


-00 

In a manner similar to the derivation of the power spectrum, the cospectrum 
can he derived as follows: 


1 * 12 ( 0 ) = / C 12 (a))do} 


00 

= lim \ C 12 ( 10 ) dw 

dco-*o ^ 


n=-oo 


00 

= lim 2 ) C 12a (gj)Aj 

Zlo->o L 
n=o 


Ris^(0) e = lim 2c 12 ti* 

AjO-* o 


C 12 a = — Rx2 a ( 0 ) 

Aj P/\m x /Ao ' ; e 

AjO-^o 


(B8) 


In general, 
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C 3 _ 2 a — lira. tTa R 12 a (O) 

. 2Aui J " :i At l r ' 


i r T 

= dim 2^ / Xl A, (t)x2 Ara (t)dt 


(B9) 


AaJ-»-0 

T^oo 


This completes the derivation for the cospectrum. Now, for the quadspectrum, 

/ OO 

Q 12 (u)) ( 


In general. 


) du) 


-Aw' 


2 Aw 



OO 

v 

iim 

dco-^o 

1 


n= -oo 


OO 
V 1 

lim 

2 y 


z_, 

n=o 

lim 

2Qi2 z 

lim 

Aa )-*0 

— i— i 

2Aw 


I dw 


; Aw' 


( 'j )Rl2 Aw^ 0) ( 


Ql2 Ao) } im 2Ao) 
ziw ->o 


(-j) R l 2Ai) (°) 


T 

= lim / Xi ^ (t)X2 ^ (t)dt 
Am — 


A 0--+0 

T OO 


i r T 

= lim / Xi A ( t ) x 2 A (t)dt 

2AwT / AjJ x ' A*r ' 
Aw-»-o °o 


T-^o 


(BIO) 


(Bll) 


This completes the derivation for the quadspectrum. 
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APPENDIX C 




DERIVATION OP OUTPUT EQUATIONS FOR ANALOG METHOD 
OF POWER AND CROSS -POWER SPECTRUM ANALYSIS 


The schematic diagram for this method is shown in figure 3. Let the 
inputs he 


00 



i=i 


osc 


B sin io 0 t 
-B cos oo 0 t 


for 

for 


0° phase 
-90° phase 


(C3) 


where 0j_, in general, does not equal cp^. The values 0^ and cp^ indicate 
that the inputs, in general, are random phase with respect to the oscillator. 
With the summation understood, the filtered outputs are 


x iAj(t) = | Ai GOS &1 ( cl <- a ) 

*i^(t) = | Aj_ sin &! (C 4 b) 

x 2 ^(t) = | Ci cos 6 2 (C 5 a) 

x|^(t) = | Ci sin 6 2 (C 5 b) 


where 


The power spectrums of 


= (o)j_t - U) Q t 
5 2 = (wj_t - u) Q t 
Xi(t) and x 2 (t) 


+ 0i) 

+ <Pi) 

in terms of mean square power are 
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* 


$;L1 ^ Aw T J 


T r 


ll Aw 


(t) + x-l^U) 


dt 


* si (Ki2 + Ks * + • 


+ K n 2) = A) . (06a) 


<J>22( w ) = 


_ 2 _ 1 
Aw T 


'T r 


x |ho (t) + x 2 w (t) 


at 


where 


= ^ (ii 2 + ^2 S + • • • In 2 ) = 


2 (Output 2) 
Aw 


A,-B 

K i “ — > 


CqB 

It = -o- > 


i = 2 , 3; 


(C 6 b) 


Wi + w 0 j > w c , |uq - w I < 


uv 


Equations (c6) have no ambiguity and are correct forms for random signals as 
well as for complex or periodic signals. The division by Aw in equa- 
tions (C 6 ) is for bandwidth normalization., and should be used for continuous 
spectra only. The cospectrum and quadspectrum in terms of mean squared power 
are 

'T r 


n l \ 2 1 

Cl2 ( w ) “ Aw T 


Xi^(t)x 2 A Jt) + x? /Vi (t)x| A , i (t) 




■Ax)' 




dt 


Ax) 


Output 3 


T r 

Qi2(u) = ~ f [ x izw (t)x 2^o (t) “ x °Aw (t)x2 2w (t) J dt 


(CTa) 


= - Output 4- (C7b) 

Ax) 

This completes the derivation for the output equations. To show that the 
quadspectrum has no ambiguity in sign as a function of frequency, 


>T ^ 


n 


-j j ' •gd. . 

Output 4 = - jp / — ) (Aq cos SiCq sin 6 2 


- Aq sin &iCq cos S 2 )dt 


(C8) 
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To simplify equation (C8), let 


Then 


xi(t) = Ai sin wit 
x 2 (t) = Ai sin(cjit - 0 2 ) 

0 2 = a variable, from 0° to 36O 0 

&1 = (till - w 0 )t 

S 2 = (wit - oj Q t - 0 2 ) 


and 


Output 4 


Ai g B 2 1 
4 T 


[cos((Ji - co 0 ) t sin(a>it - u> 0 t 


- sin(oj 1 - G) )t cos(wxt - out - 0 2 )]dt 


Case 1. GJi > by Aj q ; ui = co + Zio c 


a £t>2 j r T 

Output 4 ^ ^ / [cos(z^J Q )t sin(£w 0 t - 0 2 ) 

'-'o 

- sin(ZiJ 0 )t cos(z^o D t - 0 2 )]dt 
2-d 2 n r T 


Aj^Bf 1 

4 T 


[(-cos 2 Aj0 o t sin 0 2 ) - (sin 2 Aj 1 


2r,2 


_ A-i B 1 


sin 0 2 


Case 2. Wx = w c 


, Ai 2 B 2 1 

Output 4 = - - j^ — ^ 


.T 


( -sin 0 2 ) dt 


= ^ sin 0 2 


(C9) 

(CIO) 

- 0 2 ) 

(Cll) 


sin 0 2 ) ]dt 

(C12) 

(C13) 
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✓ P ase 3. * < W Q j COj_ = ( j) q - AJq 


A 1 P 

Output 4 - - ^ / [-cos(^w 0 )t sin(z^ Q t + 0 2 ) 

J 0 

+ sin(z^o o )t cos(Aw 0 t + 0 2 )]dt 

a 2 b 2 i r T 

^ — rji / [(-cos 2 Z^j Q t sin 0 2 ) + (-sin 2 £w Q t sin 0 2 )]dt 

Jo 

(01k) 


A l 2 ® 2 • « 

= — 1 ^ — sin 02 


Thus, for all three cases, output 4 (a quadspectrum output) is always a plus 
sine function, and there is no ambiguity. 
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